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j> ■ Abstract 

We study the chiral phase transition with staggered fermions on 
. the lattice at finite temperature in the strong coupling limit. The ther- 

\& ' modynamic potential is derived in the large d approximation where 

o" 

. the conventional method and leads to a simple physical interpretation 

for the approximation scheme. 



%*• 1 Introduction 

si ' 

The chiral phase transition has been extensively studied in the context of 
the large d expansion at strong coupling. The mean field result at finite 
temperature shows a second order phase transition 0-0. The flavor de- 
pendence and the critical exponents of the transition have also been studied 
[||. ||. However, in the conventional treatment of the large d expansion, the 
links in the temporal direction are treated quite differently from those in the 
spatial directions, leading to a complicated determinant calculation. This 
more accurate treatment of the temporal links raises a question about the 
overall consistency of this application of strong coupling method to the QCD 
phase transition. In this paper, a better motivated and simpler approach is 
presented. The final result of this new approach agrees with that of previous 
works to leading order in d. 

We start from the partition function of lattice QCD with staggered fermions. 
In section 2, we derive some basic formulae relating the chiral condensate to 



the thermodynamic potential. In section 3, we take the strong coupling limit 
g 2 — > oo. After integrating out the links in the spatial directions, we make a 
large d expansion and derive the thermodynamic potential in a saddle point 
approximation to the integration over an intermediate variable. In section 4, 
we integrate over the time links and the fermion fields to explicitly evaluate 
the thermodynamic potential. The phase transition is studied afterwards. 
Our method is much simpler than the conventional treatment, yet it gives 
the same result to leading order in the 1/d expansion. Finally, we give a 
physical interpretation of the calculation which is quite obvious in this new 
approach. 

2 Basic Formulation 

We consider lattice QCD with staggered fermions on a d + 1 dimensional 
asymmetric lattice at finite temperature. The lattice size is N s d x N t with 
N t kept fixed and N s — ► oo. The lattice spacing is a s in the spatial directions 
and at in the temporal direction. The partition function with an external 
source a(x) coupled to xx( x ) is defined as 



Z[a] = / [dUdxdx] exp{-(S G + S F ) + ^2<r(x)x(x)x(x)} , (1) 



X 



where the gauge action Sq and the fermion action Sf are 



S G 



S f 



+AE{l-^r(f/ Dt + f/t t)}) 

o f Y,ix(x)Ut{xMx + t) - x(x + t)U t (xYx(x)} 



(2) 



X 



+ Vk(x){x(x)U k (x)x(x + k)- x(x + k)U k (x) ] x(x)} . (3) 



X k=l 



In the above equations, 






(4) 
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and rjk(x) = (— l) Xl ^ hXk is the staggered phase factor. Define 

Z[a] = e w W . (5) 

The Legendre transform of the free energy W[a] is the thermodynamic po- 
tential 

rm = J2^Mx)-w[a}, (6) 

X 

where 

5cx(x) <5$(x) 
The normal chiral condensate is related to $(x) by 

(5Ma 



(x(z)x(z)) 



6 a(x] 



= Hx)Lx)=o = Hx)\srm =0 ■ (8) 

a(x)=0 S * (X) 



Thus the expectation value (xx) equals the value of $ at the stationary 
points of the thermodynamic potential 



3 Large d expansion at strong coupling 

We work in the strong coupling limit, g s — > oo, g t — > oo. To leading order 
in (3 S and (3 t) we will neglect Sq in the functional integral. We will start 
by integrating out the gauge links in the spatial directions and applying a 
self-consistent large d expansion to the fermion fields. The leading quadratic 
terms in the quark bilinear XX w iU be linearized by introducing a Gaussian 
integral over an intermediate scalar field \(x). We will then make a saddle 
point approximation to the X(x) integral and derive the thermodynamic po- 
tential in a simple form. The integrations over the time links and fermion 
fields are left to the next section. 

After integrating over all the C/^'s in the spatial directions, we get 

Z[a) = / [dUtdxdx] exp {^cr{x)x(x)x{x) 

J X 

+ ^ Y.^)U t {x) X {x + t) - x(x + t)U}(x)x(x)) 

^ x 

+ J2J2F(x,x + k)}, (9) 

x k=l 
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where 



F{x,y) = j^x(x)x(x)x(y)x(y) + 32N 2^ N _ ^ (x(x)x(x)) 2 (x(y)x(y)) 2 

+ ■■■. (10) 

F(x,y) is a polynomial of the composite fields x( x )x( x )x(y)x(y) an d 
£ iit2..*N £ jih--3N^i( x ^2( x ^ _ _ _ X* JV (^)x- ,1 (y)x j2 (2/) • • • X^ N {y) which is symmet- 
rical between the two spatially neighboring sites x and y. N is the number 
of colors. This is a result from the Uk integral and from the fact that we 
are neglecting the gauge action at strong coupling. The summation over d 
neighbors for F(x, y) in equation (§) leads to a large d expansion 0. It turns 



out that (xx) 1S on the order of J 1/d and the higher order terms in F(x, y) 
are suppressed by powers of 1/d. A naive justification for this 1/d expansion 
would be the mean field approximation. In that approximation, there would 
be a common factor d in front of all terms in F and the 1/d expansion would 
amount to a standard loopwise expansion. The actual derivation is similar 
but more involved. For our current calculation, it is sufficient to keep only 
the leading term in F(x,y). Follow the standard treatment in we define 

vM = ^T,(s a+ i^ + s x _ hf ) 1 (ii) 

which satisfies the following equations 

V(x,x') = V(x',x) , ^^(x,x')=l, Y. V ( X ^ X T 1 = 1 . (12) 

X X 

We rewrite the leading term of F(x, y) in the exponent as 

1 d _ _ 
exp{— Yl X{x)x{x)x{x + k)x(x + k)} 

4iV X k = l 

= ex P{j^ J2 x(x)x(x)V(x, x')x(x')x(x')} 

x,x' 

r N 

= j [d\] exp{- 7 £ X(x)V(x, x'r 1 ^') + Y X(x)x(x) X (x)} . (13) 
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Substituting the above equation into eq (|9|) and making the change of vari- 
ables X(x) — > X(x) — a(x), we have 



[a] = / [d\dU t dxdx\ exp X{x)x{x)xl 



X 



c 

+ ^ Y,(x{x)U t {x) X {x + t) - x(x + t)Ui(x) X (x)) 

"7 E(M^) - cr(x))V(x, x')- 1 (A(x') - a(x')) } (14) 

x,x' 

[d\] exp{ — — y^(A(x) — cr(x))V(x, x')^ 1 (X(x') — cr(x')) 

x,x' 

+ y £A(X(x))} t (15) 



where 



exp{5>(A(*))} 

X 

[dU t dxdx] exp A(s)%(a;)x(a;) 
+ I J2(x(x)U t (x) X (x + t) - x(x + t)U}(x)x(x)) } . (16) 

We now derive the thermodynamic potential. Make saddle point approxima- 
tion to the X(x) integral, 

Z[a] = e^l 

= exp {-^ Y,( X ( X ) ~ °(x))V(x, x')-\X{x') - a{x')) 

+ £^(A(*))}. (17) 

X 

The saddle point condition is 

~ E ^( x > X T (A(a/) - a{x')) + A'(A(x)) = . (18) 

(Jj 1 



We find 



$(*) = = A'(X(x)) , (19) 
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T[$] = 2<7(z)$(a:) - W[(j\ 

X 

= -^E^M^Wy) + (2°) 

And _B($(o;)) is the Legendre transform of A(X(x)) 

B($(x)) = \(x)$(x) — A(X(x)) . (21) 

4 Mean Field Results 

In this section, we integrate out the time links and the fermion fields to 
obtain A(X) explicitly. We then study the thermodynamic function r[$]. 
We choose a different integration order from the conventional treatment and 
give a simpler derivation. 

To calculate A(X), we will first integrate out the UfS, keeping only the 
leading term in the XX- Then we will introduce a scalar field Xt(x) in order 
to rewrite the XXXX term into a bilinear form. Finally we will integrate over 
the Grassmann fields XX- 

exp{5>(A(a:))} 

X 

= / [dxdx] exp {£ H x )x(x)x{x) + ttt E X(a0x(a?)x(z + *)x(^ + *) } 
= J [dxdx] exp {£ A(x)x(x)x(x) + — ^ x(x)x(x)V t (x, x')x(x')x(x) } 

x x,x' 

f N 
= J [dXtdxdx\exp{-—^2Xt(x)Vf\x,x r )X t (x') 

^ x,x' 

+ E( A (^) + *t{x))x(x)x{x)} 

X 

r N 

= J[dX t }exp{--Y<Mx)V t -\x,x')X t (x')}x]l(X(x) + X t (x)) N , (22) 

' x,x' x 

where 

V t (x,J) = l/2(5 a -, +t + 5 a -,_i). (23) 
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Making the assumption that both A(x) and Xt{x) are independent of t, we 
get 



A(X(x)) = ±- log I J dX t (x) exp{-^A,(x) 2 } x (X t (x) + X(x)) NNt j . 

(24) 

In order to study the thermodynamic function we shall find the 

Legendre transform of A(X). Let us first discuss two regions of values for A 
where we can make simple approximations: 

1) If A is large, on the order of Va, the leading term in A(X) is X NNt . So 
we have 

A(X) = N log A + const , (25) 
$ = A'(X) = j. (26) 

It is natural to set $ to a constant at this stage when we evaluate so 
that 



x,x 

d 



AN 



$ 2 + N]og$)N?N t , (27) 



(XX) = Hx)\s m=(] = NJ^. (28) 



This is the known result giving chiral symmetry breaking at zero temper- 
ature. We see that both N t and £ have disappeared from the final answer 
for (xx) j which means we are essentially working with a d + 1 dimensional 
system, identical to the original system at T = 0. The chiral condensate is 



proportional to yl/d, which justifies the large d expansion. Chiral symmetry 
is always broken. 

2) If A is small, we can expand A(X) in a Taylor series. 

A(X(x)) = -J-log{C(l + a 2 A 2 + a 4 A 4 + ...)} 

= i-{logC + a 2 A 2 + (a 4 -|)A 4 + ...}, (29) 
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where 



C = /dA t exp{-^A?}A^ 



1 ( NNt \ __ r NN t 2 NNt _ 2 _(NN t . 



«2 = ^ ^ 2 j J d\ t exp[ — — A t ]A 4 

1 { NNt \ f „ r NN t 2 NNt _t 



a4 = C\ 4 J J dX t e M — 12~\']\ 
(iViV f ) 3 (iViV f -2) 

ee 4 



(30) 



Neglecting the constant term, and taking iViVi ^> 1 just for simplicity, we 
get 



AfXf^ N " Nt ^ N 4 N? 4 

^(A(z)) = ~^^ A ~^^ A ' ( 31 ) 



£($(£)) = _J!_$2 + _^_ $ 4 (32) 

V V ;; AN 2 N t A8N*N t y 1 

Again we take to be a constant field in space, 

r[$l = [(-— + -4j— )$ 2 + — A ® 4 }NfN t . (33) 
1 J [V AN AN 2 N t J A8N 4 N t 1 s K J 

Normally, when d is large, the coefficient of the $ 2 term in the above equa- 
tion has a negative sign and we have broken chiral symmetry. When the 
temperature is raised, £ becomes large. If £ is on the order of \fd, the coeffi- 
cient can change sign and a second order phase transition occurs. Above the 
transition, (xx) = an d chiral symmetry is restored. The critical £ for the 
chiral phase transition is, to leading order in 1/d, 



a 



& = - = yjNN t d . (34) 



And the transition temperature T c is 



— ^ = -t-- (35) 
N t a t N t a s a s \ N t 
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In general, large d and high temperature are two competing factors which 
determine the phase of the system. This amounts to finding the true mini- 
mum of the effective potential in eq (|Hj), 

V e f f (X) = ^(X(x) - a(x)) 2 - A(X(x)) , (36) 

in the saddle point approximation. For a fixed d at low temperature, the 
system is in the chiral symmetry broken phase. A second order phase tran- 
sition to the chirally symmetric phase occurs when the temperature is raised 
to the order of \fd. 

Our result above agrees with the conventional calculation although our 
method is different. In the conventional treatment, A(X(x)) is calculated 
in the following order: One first integrates over the Grassmann variables 
XX, resulting in a determinant involving the links in the temporal direction. 
One then performs the integral over these UtS. Thus one performs an exact 
integration of the time links after one make the large d expansion. However, 
this final integral is quite elaborate and seems unnecessary. In this paper, we 
treat the temporal and spatial links similarly in that they are all integrated 
out at the beginning. We then do the XX integral and finally integrate over 
the intermediate field X t . We arrive at an A(X) that has the same structure 
as that given by the conventional treatment, 

A(X) = C\og(a + a 2 X 2 + --- + X NtN ). (37) 

The difference lies only in the higher order 1/d corrections to the coefficients 
in the logarithm of A(X). Also, instead of using an effective potential which 
depends on the intermediate integration field, we derive the thermodynamic 
potential which is directly related to (xx) y i a ec L ©• 

5 Discussion 

When d is large, one normally expects zero temperature physics because the 
thermal fluctuations tend to average out in this limit. In terms of QCD, this 
leads to the chiral symmetry broken phase. However, even in the large d 
limit, a chiral phase transition can occur if the temperature is high enough. 
From eq (|3~3"D , there are two terms contributing to the coefficient of the $ 2 
term in the thermodynamic potential r[$]. The first term comes from the d 
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neighboring sites which tends to drive the system into the chiral symmetry 
broken phase, and the second term comes from the thermal fluctuations in 
the Euclidean time direction. When the temperature is high enough, on the 
order of yd, the thermal fluctuations from a single time direction will be 
comparable to the total effect from all 2d spatially neighboring sites and the 
chiral phase transition occurs. In the conventional strong coupling calcula- 
tion, the role of the thermal fluctuation in the time direction is somewhat 
hidden which leads to an unnecessarily complicated derivation. In this paper, 
a simpler and more unified approach is presented. The only difference in the 
treatment of the coupling between fermion bilinears in spatial and temporal 
directions is that for the spatial direction, the intermediate integration vari- 
ables are replaced by a constant determined from a saddle point condition; 
while in the temporal direction, the corresponding quantity must be explic- 
itly integrated over to take the fluctuations into account. The chiral phase 
transition is second order in the large d limit at strong coupling. 
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